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ABSTRACT 


This thesis presents an analysis of the performance of a 
local/mobile radio communications system utiiizing the Aloha 
random access protocol with variable length packets. The 
capture phenomenon due to the near/far effect that ennances 
the performance of the system is investigated. A tagged packet 
will capture the base station if its signal-to-interference 
ratio exceeds a threshold y,. Because of the near/far effect, 
users near the base station typically have a stronger signal 
than those farther away. A multiple ring model is used to 
alleviate this problem. Users in one ring employ different 
retransmission strategy from those in other rings. A shorter 
retransmission delay is allocated to users in rings farther 
from the base station than in those closer to the base station 
in order to achieve approximately the same average delay 


throughout the network regardless of location. 
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I. INTRODUCTION 


Aloha protocol is a random access protocol in which users 
transmit their packets at will. Because of the random times at 
which users may decide to transmit, there is a chance that two 
or more users transmit at overlapping times. This results in 
a collision. Whether the overlapping covers one bit or a whole 
packet, the collided or overlapped packets will be discarded. 
Aloha protocoi has inherent characteristics of instability. As 
the traffic intensity increases, the probability of collision 
increases as well, leading to possible instabilities of the 
system. Throughput is limited to some maximum value less than 
the channe] capacity. The particular values depend on the 
access mechanism. For example, the maximum value for pure 
Aloha is 0.5e7' or 0.18. Despite its poor throughput and its 
inherent instability, the Aloha protocol is still a prime 
candidate for many local and mobile packet radio systems 
because of its ease in implementation. In many practical 
applications, Aloha performance in local and = wmobile 
communication systems is affected by the near/far effect as 
well as multi-path fading. The near/far effect arises when 
users are located at different distances from the base 
station. The power levels of the received packets may vary 
substantially due to spatial attenuation of the signals, and 


the packet with highest power level has the best chance of 





being received correctly. Multi-path fading effects also cause 
random power variation among the received packets. Hence, th? 
near/far effect and multi-path fading combine to allow a 
receiver to correctly receive a packet even when multiple 
packets are transmitted simultaneously. Recent research has 
shown that each of the above effects do improve the throughput 
of an Aloha channel considerably [1-3]. However, this is 
achieved in favor of the users near the base station and at 
the expense of the far users. The far users have to retransmit 
their packets more often than the near users and hence suffer 
a longer packet delay. 

To correct this problem, a ring model in which the 
coverage area is divided into 3 concentric rings with the base 
station in the center is investigated in this thesis. Users in 
different rings use different strategies to retransmit their 
collided packets. A ring that is far from the base station 
will have a shorter average retransmission delay than a ring 
that is closer to the base station to compensate for the 
higher capture probability of the near users. As a result, 
users in a ring model may achieve approximately the same 
packet delay. 

This thesis presents an analysis of the near/far effect on 
channels employing variable packet Aloha in a local/mobile 
radio ring model. Unlike idealized Aloha in which a tagged 
packet is only received successfully when it does not collide 


with other packets, the capture effect enables the base 








station to accurately decode the tagged packet even when it is 
involved ina collision. The result of the analysis is used to 
find the equal riny-throughput of a network consisting of 
three concentric rings. 

Chapter II presents a typical system model of local/mobil 
radio communications. Chapter III analyzes the dynamics of 
number of interfering packets in variable packet Aloha 
channels. The probability of capture of a tagged packet by the 
base station and the system throughput will be established in 
chapter III. Finally, in chapter IV a network of three rings 
is analyzed to find the ring widths such that the throughput 


in each ring is approximate equal. 














II. SYSTEM MODEL 


The system model is a network consisting of users assumec 
to be uniformly distributed in an annular ring surrounding che 
base station. Figure 1 depicts the system model. Assume ‘hat 
the user traffic density is y(:), where r is the distance to 
the base station. The density g(r) can be thought of as the 
normalized packet traffic per unit area at a distance r from 
the base station. The pavket generation rate in a given area 
depends on r and is independent of direction. Then the traffic 
in a differential ring of width dr at a radius r from the base 
station is [4] 


G(dr) =2nrg(r) dr (1) 


The traffic generated within a distance r of the base statior 


is 


G(r) = an {xg(x) dx (2) 
0 


The total offered traffic in the network can be obtained 


as G = G(r ~ ©) 


G= an {xg(x) dx 
0 








The spatial distribution of packet generation is found to 
be the probability a packet is generated within a distance 


Ree 


PG 
F,(r) = Pr{Psr} = 2h [xg (x) dx (4) 

0 
Ther fore, the probability density function for packet 


generation is found to be 


f(r) = 42 rg(r) (5) 


z 





Figure 1. The system model 


Assume that users are uniformly distributed in an annular 


ring R, < r = R, around the base station, then 


G 
—————_ , 2, s7s 
g(r) | n (R2-R2) 2 R, (6) 
0 , otherwise 


Substituting (6) into (5), we obtain 


2r 
fp(r) =) RP - R 
0 , otherwise 


»R srs k, (7) 











Figure 2 shows the distribution of packet generation of a 


ring with R, <r < R,. 


g(r) 


n (RZ - RZ) 


R, R, r 


Figure 2. User traffic density g(r) 


There are two main effects that influence the performance 
of local/mobile radio communications, namely the near/far 
effect and multi-path fading. To model the near/far effect, 
assume that signal power depends on the spatial distribution 
of users in the network. The mean power of a packet at a 


distance r from the base station is of the general form 


w= cr? (8) 


where a is between 3 and 5 for ultra high frequency (UHF) 
applications. The constant c depends on transmitted power and 
the gains and the heights above the ground of the base station 
and mobile antennas. Assuming all users are identical we can 
normalize the constant c to unity for all users. Hence, 


were (9) 


Let the power random variable W be defined by 





W= R-¢ (10) 
The event {W < w) occurs when event A = ({R™% < w) occurs. 
A can be rewritten as A = {R > w'/*"). Thus, 
Py(w) = P[R2w/*] (11) 
Therefore, (11) is equivalent to 
F,(w) = 1- F,(wii/*) (12) 
Differentiate both side of (12) 


-1/a 
f,(w) = ~ £y(wov/e) GOW) (13) 


Expanding (13), we get 


md 


—-l 


ae r 
£,, (w) = fy (w ©) 2 (14) 


Substituting (7) into (14), we obtain 


-(241) 
f,(w) = —__—4__w : (15) 
a (R3 = Ry) 
For a typical value a@ = 4 in UHF communications, (15) 
becomes 

- (3/2) 
f,(w) = —____ (16) 

2(Rz - Ry) 


One of the important factors in determining the quality of 
signal reception is the number of interferers to the tagged 
packet. Thus, it is of interest to analyze the dynamics of the 


system. This is presented in the next chapter. 





II. SYSTEM DYNAMICS 


In this study, we propose variable packet Aloha instead of 
the conventional fixed length packet Aloha. Variable packets 
do not require padded bits to fill up a packet to its 
specified length which may reduce the throughput of fixed 
packet Aloha. Variable packets also do not require bit 
stuffing and de-stuffing logic which may increase the cost of 
hardware. The channel traffic rate is assumed to obey a 
Poisson process with parameter g packets per second. The newly 
created packets are of variable length x (seconds) and are 
independently and identically distributed according to the 
probability density function a(x) with mean x. All 
retransmitted packets have their length x redrawn afresh from 
a(x) for every retransmission. 

The probability that n packets interfere (overlap) with an 
arbitrary tagged packet during its transmission is given by 


[5] 


Pr{n} = AgGeul’ e-g(x+u) (17) 


where x is the tagged packet length, u is the random length of 
a preceding packet, and n is the number of interfering 
packets. 

The probability that the tagged packet encounters no 


overlap from previous packets and the next inter-arrival is 








larger than x is given by n = 0. This is also the probability 


of successful tagged packet transmission for the case where no 
interfering packets can be tolerated by a tagged packet. 
This analysis considers the case when the base station is 
captured by a tagged packet in the presence of interfering 
packets if the ratio of the tagged packet signal power and the 
joint interfering packets signal power exceeds a threshold Yo° 
The number of interfering packets in a variable packet 
Aloha channel is a stochastic process which is driven by four 
types of interferers. The first type of interferers are 
designated as Early-Late packets. These are packets that begin 
before the tagged packet and end after the tagged packet. The 
remaining types of interferers have similar descriptive 
nomenclatures and are designated as Early-Early packets, Late- 
Early packets, and finally Late-Late packets. Figure 3 depicts 


four types of interferers relative to the tagged packet. 


Figure 3. Four types of interfering packets 








The arrival and departure of n interfering packets 
partition the tagged packet into several non-overlapping 
intervals of random length. The number of these intervals, 
given a known number of interfering packets n and their 
corresponding type of interfering packets found in is [6] 


Total Intervals = number of Early-Early packets 
+ 2 x (number of Late-Early packets) (18) 
+ number of Late-Late packets +1 


The relationship in (18) is based on the following facts. 
There is initially one interval, the tagged packet itself, 
when there are no interfering packets. The presence of a 
Single Early-Early packet divides the tagged packet into two 
smaller intervals; it adds one more to the total number of 
intervals. The presence of a single Late-Late packet also 
partitions the tagged packet and provides one more to the 
total number of intervals. The presence of a single Late-Early 
packet will further divide the tagged packet and adds two more 
to the total number of intervals. The Early-Late packet 
overlaps entire duration of the tagged packet and therefore 
contributes no extra intervals to the total number of 
intervals. 

Each interval has associated with it a total number of 
interferers I(t,), where t, is the start time of the i'th 
interval of the tagged packet. Each unique sequence of total 
interferers per partition I(t,;) for all i is called a 


realization. Figure 4 depicts a specific realization. 


10 








Figure 4. An example of a specific realization. 
The total number of interferers I(t,;) = (2,3,4,3,2}) 


For this specific realization, the number of interfering 
packets is four. There is one of each type of packet: one 
Early-Late packet, one Early-Early packet, one Late-Early 
packet, and one lLate-Late packet. The total number of 
intervals during the tagged packet according to (18) is five. 
The total number of interferers for each interval is shown at 
the top of Figure 4. We are interested in finding every 
possible realization for a given n users. From those 
realizations, the number of realizations C; in which the 
maximum number of interferers present at any time during the 
tagged packet interval does not exceed some number j, will be 
calculated. 

Finding all possible realizations is very tedious but is 
easily obtained via computer. TABLE I lists the algorithm of 


the main program. 
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TABLE I. MAIN PROGRAM ALGORITHM 


% N+1l is a number of users in the network 
for n = 1:N, 
% call subroutine make_comb to find all possible 
% combination of interfering packets. 
length = make_comb() ; 
for i = 1:length, 
call subroutine realization to build a binary 
tree of all possible realizations 
realization(); 
call get_cj to count number of realizations 
which has maximum number of interferer not 
exceeding j 
get_cj(); 
end 
end 


0 dP dO de IP 


For a given n from 1 to N, where N is the maximum number 
of interferers in the system, subroutine make_comb finds all 
possible combinations of four types of packets taken n at a 
time. Basically, subroutine make_comb assigns different values 
each time for each type of interfering packets such that the 
total of all interfering packets is equal to n and enumerates 
the number of distinctive combinations found. Subroutine 
make_comb uses the algorithm found in [7] to generate all the 
combinations. For example, given n = 2, the set of 
combinations of four types of interfering packets will be U = 
{(2,0,0,0), (0,2,0,0), (0,0,2,0), (0,0,0,2), (1,1,0,0), 
(1,0,1,0), (1,0,0,1), (0,1,1,0), (0,1,0,1), (0,0,1,1)}. 

For each combination of four packet types taken n at a 
time, the subroutine realization is called to find the total 
number of realizations. The root of the binary tree 


corresponds to the first interval and contains the initial 
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number of interferers. In the first interval, the interferers 


can only be those packet types beginning before the tagged 
packet; i.e., Early-Late packet(s) and Early-Early packet(s). 
The initial number of interferers will be the sum of those 
packets. The children of the root of the binary tree 
corresponds to interval 2 and their children corresponds to 
interval 3 and so on for the rest of all intervals. The next 
intervals of the tagged packet contain either one packet more 
or one packet less than the number of interfering packets of 
the previous interval because their existence depends on 
either the arrival or the departure of one of the interfering 
packets. However, the last interval should be equal to number 
of those packets ending after the tagged packet ends. We 
choose to implement the binary tree by having the left child 
node taken care of the case of arrival of an interfering 
packet and the right child node for the case of departure of 
an interfering packet. In other words, the number of 
interferers stored at a left child node of a parent node is 
equal to the parent's number of interferers plus one if the 
parent node is not greater than n, the maximum number of 
interferers. If it is greater thann , that left child node is 
set to NULL; i.e., any realization containing that node is 
invalid. The number of interferers stored at a right child is 
equal to the parent node's minus one if the number of 
interferers of the parent node is greater than the minimum 


number of interferers which is equal to number of Early-Late 
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packet(s). If it is not greater than the minimum number of 
interferers then the right child node is set to NULL. Any 
realization which has the number of interferers in its last 
interval not equal to the sum of the number of Early-Late 
packets and Late-Late packets is also invalid. Figure 5 shows 
a binary tree of a specific combination of four types of 


interfering packets (1,1,1,1). 





OOOW WO @ © W 


Figure 5. A binary tree implementation of (1,1,1,1) 


The subroutine get_cj traverses the binary tree to collect 
all the valid realizations. The count of realizations which 
have the same number of maximum number of interferers j at any 
time is C;- 

TABLE II shows the output of the computer program for n =1 
to 7. For the case of n = 1, there are only four possible 
realizations, namely the interferer can be Early~Early, Early- 


Late, Late-Early, or Late-Late packet. Each realization has 
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only one known interferers. Therefore, the value of C,(1) is 
four. Other values of C;(2) such as C,(1), C,;(1), etc. are 0. 
For the case of n = 2, there are 10 combinations of four 
packet types taken two at a time. These combinations generate 
14 unique realizations. Of these 14 realizations, four have a 
maximum number of interferers of 1 at any time, and 10 have a 
maximum number of interferers of two at any time. The value of 
C,(2) is 4 and C,(2) is 10. Other C; (2) values are 0. The 
function I, is the total number of all possible realizations 
for n known interfering packets. The function f(n) is the 


inverse of I,. 


TABLE II. MAXIMUM NUMBER OF INTERFERING PACKETS GIVEN n 


100 
200 
400 1698 2250 1565 
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III. PROBABILITY OF CAPTURE AND THROUGHPUT OF VARIABLE PACKET 


ALOHA CHANNELS 


A. PROBABILITY OF CAPTURE 

As we have mentioned earlier, the near/far effect arises 
when users are scattered at different distances from the base 
station and causes the power levels of received packets to 
vary substantially. Thus, the arriving packet with the highest 
power level has a best chance of being received correctly. 
This model in which a variable packet Aloha protocol is 
employed has better throughput than the ideal Aloha protocol 
due to the near/far effect. A tagged packet has the capability 
to capture the base station even when it is involved ina 
collision. We assume that a tagged packet will capture the 
base station if its signal-to interference ratio exceeds a 
threshold y,. We also assume that thermal noise is negligible 
since it is small compared to the competing interference of 
interfering packets. In this section we calculate the 
probability of capture for the variable packet Aloha ring 
model and develop an expression for the throughput. 

Let the random variable X denote the power of a tagged 
packet with the probability density function (pdf) f,(x), and 
the random variable Y denote the joint interfering power of n 


> 1 interfering packets with the conditional probability 
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density function f,(y|n). With the assumption of independence, 


the conditional pdf f,(y|n) is found by the relationship [6] 
7 2 
fy(y|n) = £(n) YC, (n) [£, (w) 17 (19) 
jel 


where j@ denotes j-fold convolution, f(n) is the inverse of 
total number of realizations given n interfering packets, and 
C;(n) is the number of realizations in which a maximum number 
of interfering packets equals j givenn interferers. Both f(n) 
and C;(n) are listed in TABLE II. 


Let the random variable Z be the signal-to-interference 


ratio 
Z= 4 (20) 
then the probability density function of Z is [8] 
f,(z2|n) = [vtuy (x yln ay (21) 
Q 


Given the fact that X and Y are independent random 


variables, we obtain 


f,(z|n) = fy tx(v2 fy (yla) dy (22) 
0 


The probability of capture, that is the probability that 
Z is greater than a threshold y, given n interfering packets 


Pr{capture|n} = Pr{y, < Z < ~|n} (23) 


is obtained by integrating f,(z|n) from y, to ©. Hence 
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Pr{capture|n} = ff, (2|n) dz (24) 
Ye 
The overall probability of capture of a tagged packet is 
found by multiplying the conditional probability of capture 
given n interfering packets by the probability that n packets 
interfere with a tagged packet during its transmission, and 
summing over all n. 


N 
Pricaptie} = > Pr{capture|n} Pr{n} (2 
n=#0 


in 
se 


where N+1l is the number of users in the system. 


Substituting (17) into (25), we obtain 
= [g(x+u)]? 
Pr{capture} = > De sine) Pr{capture|n} (26) 
n=0 % 


For the purpose of illustration, consider a network with 


R, = 0 and R, = 1. The pdf of packet power f, in (16) becomes 


fy(w) = Sw 7a) (27) 


Equation (21) and the - nditional probability of capture 
F -{capture|n} with the assumption of Y, being equal to 6 dB is 
obtained numerically using Pro-Matlab (see Appendix A for the 
program). TABLE III shows the conditional probability of 
capture for n from 0 to 7. 

In the next section we will incorporate results in TABLE 
III to find the throaghput of a variable packet Aloha channel 


in a ring model. 
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TABLE III. THE CONDITIONAL PROBABILITY OF CAPrTUhE 


Pr{capture|n} 


1.0000 


0.2002 


0.1503 


0.1202 


0.1018 


0.0894 


0.0806 


0.0740 





B. THRCUGHPUT OF A VARIABLE PACKET ALOHA RING MODEL 
The channel throughput S, in packets per mean packet time 


1S Ss 


Sc: 


ote, 8 


[Hox u) a(x) a(u) dxcu (28) 
is) 


where a(x) and a(u) are exponential density functions with 
mean x of the tagged packet length and the preceding packet 
length, respectively. The quantity ¢(x,u)/2 is the conditional 


channel departure rate [5]. 


+6(x,u) = $9(x+u) Pr(capture) (29) 


Substituting (26) and (29) into (28), we have 


19 











ae 


foatx+u) > e ~ateeu) Lg esta) T° (30) 
o n=O n 


. Pr{capture|n} a(x) a(u) dxdu 
Utilizing the binomial identity 


n 


(x+u)? = yen (322) 


k=0 


and expanding (30), we obtain 


ee 
= gxe-* eT + gue"M e-& | I” pricapture|n} 

ae pe a (32) 
23 (i)** u®-* a(x) a(u) dxdu 


Further manipulating terms in the above equation gives 


~ i g” 
=o par: 2 pr{capture|n} 


=0 

= n k+ - r -Kkp- 

lel * le-® a(x) dx. lace e-a(u) du (33) 
= n\f n-k+ - r ka- 

“ele te-Ma(u) du. [x e-% a(x) dx 


=0 


tay 


We assume that the packet length is exponentially 


distributed and the average length is x. Then we have 


a(x) = te-*/* (34) 
x 

a(u) = +e-w/* (35) 
x 


Define the average traffic rate G as 


G = gx (36) 


We apply (34), (35), (36) and [9]. 








ae e- ¥* dx = nip (st) (37) 
0 


to solve (33) 


The first integral term is evaluated as 


7 ~(g+H)x 
[ x06) e-ma (x) dx = =f xe we dx 
ie) 0 


2 (k+1) ! Ges (38) 
(ke) x 

- +1)! ~ (k+2) 

"ee (kel) (G+1) : 





Similarly, the second integral is found to be 


= 2 a -(g+)u 
fae *) @ Ma(u)du = =f ui Ke t'"du 
0 0 


a ! -(n-k+1) 
= (n-k)! (g++) (39) 
ne xX 
' (aa) | -(n-k+1) 
"Say 
The third integral term is 
futren ew a(u) du = (n-kK+1)! (G44) -(a-ke2) (40) 
s x7 (a-k+1) 
The final integral term is reduced to 
[ xre-a(x) dx = EL (G41) ten) (41) 
x 
0 


Combining all integral terms in (38), (39), (40), and (41) 


into (33), we get 


a ae n (G+1) (9+) ut 
S = = os  Pricapture|n} “a — n!(n+2) (42) 
neo 844 k=0 
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Applying (36) and the equality 


Yo n! (n+2) = nt (n+1) (n+2) (43) 
kx 


and further simplifying (42), we obtain 


oo > (n+1) (n+2) 


= Pr{capture|n}G"**(G+1)~(*}) (44) 
n=0 


For the network with R, = 0, R, = 1, 


and the conditional 
probability of capture in Table III, throughput is shown in 


Figure 6. 
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Figure 6. Throughput of a ring model 
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IV. APPROXIMATE EQUAL RING THROUGHPUT 


Tm this section, we expana cur analysis to a system of 
multiple rings. As we mentioned before, the probability of 
capture favors near users at the expense of far users 
resulting in lower throughput for the far users. A multiple 
ring system compensates for this deficiency by using separate 
retransmission strategies for users in different rings. A user 
in a ring that is far from the base station is given a shorter 
average transmission delay than a user ina ring closer to the 
base station to compensate for a higher capture probability of 


the near ring. 


A. PROCEDURE TO FIND EQUAL RING THROUGHPUT 

A procedure analogous to the procedure found in [3] is 
used to find the ring width for each ring such that the 
throughput in each ring is approximately equal. The procedure 
consists of the following iterative process: 

1. Select initial ring widths R; = {r; < r < 1r;,,,}, where 
SS Dy 23 hg Ms 

2% Find the conditional ring capture probability 
Pr(capture|n,R;}(see Appendix B). 

3. Find the traffic rate G, in each ring 


G,=GIF,(r;,,) -Fe(r;)] (45) 
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where i = 1,2,...,M and F,(r,;) is given in (B3) 


4. Find the ring throughput fer each ring from (Appendix 


Cc). 
S,= Y i. 2) (22) pr {captuce|n, Ry} 6,6" (G+1) ee (46) 


5. If S, = S, ¥...% S, for 0 ¢ S; < min{max S;} go to next 
step. Otherwise, go back to step 1. For our computational 
purpose, the maximum relative difference of +0.05 is 
acceptable for the approximate equality. 

6. Calculate the total throughput of the model S = iS;. 

For illustrative purpose, consider a system of three rings 
R, = {0 < r < 0.52}, R, = {0.52 < r < 0.8} and R, = {0.8 < res 
1}. The threshold y, is taken to be 6 dB. Fig 7 shows the ring 
throughput of ring R,, R,, R; and the total ring throughput of 
the model. The ring throughputs are approximately equa] for 0 


< Gs 0.25. 


B. AVERAGE RETRANSMISSION DELAY 

To determine the average retransmission delay T, for ring 
R,;, we need first to define a retransmission strategy when a 
collision occurs. One suggestion [10] is to choose an 
arbitrary time interval and select a uniformly distributed 
random retransmission time within that interval. Let the time 
interval covers K; message-unit times where each message-unit 


is an average x units of time. Then retransmission takes place 


within 1 to K,; time units after a collision has occurred. The 
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Figure 7. Approximate Equal Ring Throughput 
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approximate number of retransmissions is found to be 
oe 1-Pr{capture|R;} (47) 
Pr{capture|R;} 
where Pr({capture|R,;} is the probability of capture in ring R, 


per mean packet time is found by (Appendix D) 


Pr{capture|R,} = }> Pr{capture|n, R,} (n+1)G"(G+1)~(**) (48) 


n=0 
The average retransmission time, ignoring propagation 
delay, is approximately 


T,; = 0.5E(K,+1) (49) 


By fixing the value K, of the outermost ring, the 
remaining value K,; can be obtained by letting T,; = T, and 
substituting (47) into T,; and Ty expressions, we obtain 


Pr{capture|R,} (1-Pr{capture|Ry}) (Ky+1) _ 


K,; = 50 
2 Pr{capture|R,} (1 -Pr{capture|R,} ) - ee 


TABLE IV shows the number of retransmission K, and K, of 


ring R, and ring R,, given K, = 8 of X time-units. For 0 <G < 


Giax = 0-25, the average K, is 29.37 and K, is 8.70. 











TABLE IV. VALUES OF K; 


27.9389 8.6379 
28.6474 8.6694 
29.3634 8.7003 
30.0863 8.7306 
30.8158 8.7603 








Vv. CONCLUSION 


We have evaluated the probability of capture taking into 
account the near/far effect. These results are then used to 
compute the throughput of a ring model utilizing the variable 
Aloha protocol in a local/mobil radio communications. We have 
investigated a ring model in which different retransmission 
strategies are used to alleviate the near/far effect problem. 
By assigning shorter average retransmission delay to users 
farther from the base station and longer average transmission 
delay to users nearer to the base station, all users may 
achieve approximately the same average retransmission delay. 
In this thesis, we only investigated a system of three rings. 
Further work is needed to determine the number of rings that 
produces the maximum throughput. The procedure for finding 
equal ring throughput is a process of trial and error and 
requires greater effort when a number of rings is greater than 
three. The capture measure using signal-to-interference 
threshold provides a less accurate prediction of channel 
throughput than does a capture measure using packet correction 
reception probability [1]. However, employing packet 
correction reception probability as a capture measure for 
variable packets is very difficult. It is an unwielding task 
to find the probability of bit error for variable length 


packets. 
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A. 





APPENDIX A. COMPUTER PROGRAM LISTINGS 


PROGRAM TO FIND THE MAXIMUM NUMBER OF INTERFERING PACKETS 


The program is written in Unix C to find the number */ 
of realizations Cj which has the maximum number of */ 
interferers equal to j. It first calls subroutine */ 
make_comb to find all the combinations of four types */ 
of interfering packets taken n at a time. Then, for */ 
each combination it calls subroutine realization to */ 
build a binary tree of realizations. */ 

Finally, subroutine get_cj is called to find Cj. */ 


#include <stdio.h> 

#define comb_lgth 230 

#define type 4 

#define cj_lgth 10 

#define INVALID -1 

#define NODE_SIZE sizeof(struct node) 


/* 


Global variables declaration */ 


struct node 
{ 
int num_intfr; /* number of interfering 
packets */ 
struct node *left; /* left pointer of a node ina 
binary tree */ 
struct node *right; /* right pointer of a node in 
a binary tree*/ 
i 


FILE *file; /7* file descriptor */ 
int num_user; /* number of users in system */ 
int init_intfr; /* number of intefering packets in 


the first interval i.e. number of 
early beginning packets */ 
int final_intfr; /* number of late ending packets */ 


int max_intfr; /* number of maximum possible 
interfering packets i.e. number of 
users */ 

int min_intfr; 7* minimum number of interfering 
packets */ 


/*i.e. number of early-late packets */ 
int num_interval; /* number of intervals of tagged 
packet */ 
int cj[(cj_lgth][cj_lgth); /* array of Cj */ 
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7* Function declaration */ 

int make_comb() ; 

struct node *realization(); 
void get_cj(); 

void print_out(); 

char *malloc(); 


main() 


ant i,4: 

int length; /* length of array comb */ 

int comb[comb_lgth][type);/*array of combination of 
possible packet types*/ 

struct node *root; 


file = fopen("datafile","a"); 

root = (struct node *)malloc(NODE_SIZE); 
root~->left = NULL; 

root~->right = NULL; 

/* clear array cj */ 

for (i=0;i<cj_lgth;i++) 


for (j=0;j<cj_lgth;j++) 
cjJ[1}(3] = 0; 


} 
for (num_user=1;num_user<=10;num_user+t+) 


{ 
length = make_comb(comb) 7 
for(i=0; i<=length; i++) 
{ 
init_intfr = comb[i][0]+comb[i][1]; 
final_intfr = comb[i][0]+comb[i][type-1]; 
max_intfr = num_user; 
min_intfr = comb[i][0]; 
num_interval = comb[i][1] + 
2*comb[i][{2]+comb[i][type-1]+1; 
root->num_intfr = init_intfr; 
root = realization(root) ; 
get_cj (root); 
) 


) 
print_out(); 
fclose(file) ; 

} 


/* The subroutine make_comb is to find all possible */ 
/* combinations of four types of interfering packets */ 
/* taken n at a time. */ 


int make_comb (comb) 
int comb[comb_lgth] [type]; 
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int pkt_tvpe[{type]; 
int length = 0; 

int index_track = ~1; 
int num_track; 

int i, 33 


pkt_type[0] = num_user; 
num_track = num_user; 
comb(0J}[0] = pkt_type[0]; /* first combination */ 


for (i=l; i<type; i++) 
pkt_type[i] = comb[length][i] = 0; 
while(pkt_type[type-1] != num_user) 
{ 


if(num_track >1) 
index_track = -1; 
index_track++; 
num_track = pkt_type{index_track]; 
pkt_type[index_track] = 0; 
pkt_type[0] = num_track - 1; 
pkt_type(index_t.rack + 1]++; 
length++; 
for(j}=0; j<type; j++) 
comb(length)(j] = pkt_type([j]; 


) 
return(length) ; 
} 


/* Subroutine realization employs subroutine build_tree */ 
/* to build a binary tree of realizations. */ 


struct node *realization(ptr) 
struct node *ptr; 
{ 
int interval; 
int tree_level; 
struct node *build tree(); 
void build_node(); 


interval = 1; 
tree_level = 1; 
while (interval<num_interval) 


{ 
ptr = build_tree(ptr,interval,tree_ level); 
interval++; 


} 
return(ptr); 


) 
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/~ Subroutine build_tree searches for the valid paths of 
realizations */ 


struct node *bu’ld_ tree(ptr,interval,tree_level) 
struct node *ptr ; 
int interval; 
int tree_level; 
{ 
if (tree_level<interval) 


{ 
if((ptr->left)->num_intfr != INVALID) 
{ 
tree_level++; 
© r= > 1 e f t = 
build _tree/ptr->left, interval,tree_level) ; 
tree_level--; 


} 
if ((ptr->right)->num_intfr != INVAL1D) 
{ 


tree_level++; 
ee = > SF Lf Gg Hw = 
build_tree(ptr->right,interval,tree level); 
tree_level--; 
} 
} 


else 
build_node(ptr,interval) ; 
return(ptr) ; 
} 


/* Subroutine build node is to fill nodes of a binary */ 
/* tree with the number of interferers. A left child */ 
/* node has the number of interferers equa to the */ 

/* parent's plus cne if the parent's is not greater */ 
/* than n, the maximum number of allowed interferers. */ 
/* If the parent node is greater than n, the left */ 

/* child node is set to NULL, i.e. invalid. A right */ 
/* child node will be equal to the parent's number of */ 
/* interferers minus one if the parent's is greater */ 
/* than the minimum number of interferers, which is */ 
/* equal to number of Early-Late packets. */ 

/* If it is not, the right child is */ 

/* set to NULL. The final interval should have the */ 

/* number of interferes equal to the sum of number */ 

/* of Early _Late and Late-Late packets. */ 

/* If it is not, that realization is invalid. */ 


void build_n vde(ptr,interval, 
struct node *ptr; 
int interval; 
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if(ptr->left == NULL) 


{ 
if((ptr->left = (struct node *) malloc (NODE_SIZE) ) 
== NULL) 


{ 
printf ("ERRORI!!!"); 
exit(1); 
} 
(ptr->left)->left = NULL; 
(ptr->left)->right = NULL; 


} 
if (ptr->right == NULL) 


{ 
if ((ptr->right = (struct node *)malloc(NODE_SIZE) ) ==0) 


{ 
printf ("ERROR2!!!"); 
exit(1); 


} 
(ptr->right)->left = NULL; 


(ptr->right) ->right = NULL; 


} 
/* install left child node */ 
if (ptr->num_intfr<max_intfr) 


{ 
(ptr->left)->num_intfr = ptr->num_intfr + 1; 


/* invalid final interferers */ 
if ((interval==num_interval-1) &&((ptr->left) ->num_intfr 


!=final_intfr) ) 
(ptr->left)->num_intfr = INVALID; 


} 


else 
(ptr->left)->num_intfr = INVALID; 


/* install right child node */ 
if(ptr->num_intfr>min_intfr) 
{ 
(ptr~>right)->num_intfr = ptr->num_intfr - 1; 
/7* invalid final interferers */ 
if ((interval==num_interval-1) &&((ptr->right) ->num_intfr! 


=final_intfr) ) 
(ptr->right) ->num_intfr = INVALID; 


} 


else 
(ptr->right) ->num_intfr = INVALID; 
} 


/* Subroutine get_cj is to count the number of */ 
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/* realizations which has the maximum number of */ 
/* interferes not greater than j. */ 


void get_cj (ptr) 
struct node *ptr; 
{ 
int tree_level; 
int max; 
void count_cj() 


me 


tree_level = 1; 

max = ptr->num_intfr; 

count_cj(tree_level,max,ptr) ; 
) 


void count_cj(tree_level,max,ptr) 
int tree_level; 

int max; 

struct node *ptr; 


int old_max; 
if (tree_level<num_interval) 


{ 
if((ptr->left)->num_intfr != INVALID) 

{ 
tree_level++; 
old_max = max; 
if ((ptr->left) ->num_intfr>max) 

max = (ptr->left)->num_intfr; 

count_cj(tree_level,max,ptr->left) ; 
tree_level--; 
max = old_max; 


} 
if((ptr->right)->num_intfr != INVALID) 
{ 
tree_level++; 
count_cj(tree_level,max,ptr->right) ; 
tree_level--; 
} 
} 
else 
cj (max_intfr-1] [max-1]++; 
return; 


} 


void print_out() 


int total_rlzt; 
double fn; 
int: i,j? 

for (i=0;i<cj_lgth;i++) 





fprintf(file,"%td\n",i); 
total_rlzt = 0; 
for (j=073<cj_lgth; j++) 


fprintf(file,"ta " ,cj{ij(j]): 
total_rlzt += cj[{i)(j): 
} 
fn = 1.0/ (double) total_rizt; 
fprintf (file, "%d %8e\n",total_rlzt,fn); 
} 
} 


B. PRO-MATLAB M-FILE FOR THROUGHPUT OF A RING MODEL 
The following M-file is written to run on Pro-Matlab to 


find the probability of capture and throughput of a ring model 


with R, = 0 and R, = 1. 





% set w from 1 to infinitive 

w = 1:0.1:52.2; 

resolution = 0.1; 

range = 1/2; 

% find the pdf of power of a packet 

fw = range.*(w.%*(-3/2)); 

% find j-fold convolution of interfering packets 
CONV(1,:) = fw; 

CONVFW = fft(fw); 

for i= 2:10, 

CONV(i,:) = resolution *ifft(fft(CONV(i-1,:)) .* CONVFW); 
end 

% get f(n) and Cj from data file 

data 


tmp = tmp + (cj(i,k) .* CONV(k,:)); 


fy(i,:) = fn(i).* tmp; 

end 

% find fz in fz(z/n) = f2*2z*(-3/2) 

for i = 1:10, 

sum = range * (w(1)*(-0.5)) * fy(i,1); 

for j = 2:2:512, 

sum = sum + 4*range*(w(j)%*(-0.5))*fy(i,j); 
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end 


for j = 3:2:511, 

sum = sum + 2*range*(w(j)*(-0.5))*fy(i,j); 
end 
sum = sum + range*(w(513)%*(-0.5))*fy(i,512): 
fz(i) = (resolution/3) *sum; 
end 


% find the probability of a tagged packet 
pr = 2 * (6%(-1/2)) .* £2; 


% Throughput 

g = 0:0.1:5; 

for i = 1:51, 

s(i) = g(i)/((g(i)+1)*%3); 

for n = 1:10, 

s(i) = s(i) + (n+1)*(n+2)/2*pr(n) *(g(i)*(n+1)) 
* ((g(i)+1)*(-(n+3)))3 
end 
end 


C. PRO-MATLAB M-FILE TO FIND EQUAL RING THROUGHPUT 
This M-file is written for Pro-Matlab to find the equal 
ring throughputs in chapter 4, 


% set w from 1 to infinitive 
w= 1:0.05:26.6; 
resolution = 0.05; 
range(4) = 1/2; 
% pdf of power of an interfering packet 
fw = 1/2.*(w.*(-3/2)); 
% calculate the j-fold convolution of interfering 
% packets 
CONV(1,:) = fw; 
CONVFW = fft(fw4); 
for i= 2:10, 
CONV4(i,:) = resolution.* ifft(fft(CONV4(i-1,:)) .* 
CONVFW) ; 
end 
% get f(n) and Cj from data file 
data 
% find fy 
for i = 1:10, 

for j = 1:513, 

tmp(j) = 0; 

end 
for k = 1:i, 

tmp = tmp + (cj(i,k) .* CONV4(k,:)); 
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end 
fy(i,:) = fn(i).* tmp; 
end 
% polynomial fit for fy(n) 
for i= 1:10, 
fit(i,:) = polyfit(w,fy(i,:),9); 
end 
% set range of ring 1 from 0 to 0.52 
r1(1) = 0; 
r2(1) = 0.52; 
range(1) = 1/(2*(0.52%2)) ; 
% set range of ring 2 from 0.52 to 0.8 
r1(2) = 0.52; 
r2(2) = 0.8; 
range(2) = 1/(2*(0.8*%2 -0.52%2)); 
% set range of ring 3 from 0.8 to l 
r1(3) = 0.8; 
r2(3) = 13 
range(3) = 1/(2*(1 - 0.8%2)); 
% find fz in fz(z/n) = £2*z*(-3/2) 
for ring = 1:3, 
for i = 1:10, 
for z = 6:.1:24.4, 
if ring == 3, 
start = 1; 
else 
start = r2(ring)*(-4)/z; 
end 
if ring ==1, 
ends = 26.6; 
else 
ends r1(ring)*(-4)/2; 
end 
temp = 0; 
if start <1, 
start = 1; 
end 
if ends > 26.6, 
ends = 26.6; 
end 
if ends > 1, 
step = (ends - start)/8; 
y = start:step:ends; 
clear f; 
f = polyval(fit(i,:),y)? 
sum = range(ring) * (y(1)*(-1/2)) * £(1); 
for j = 2:2:8, 
sum sum + 4*range(ring) *(y(j)*(-0.5)) *£(j): 
end 
for j = 3:2:7, 
sum sum + 2*range(ring) *(y(j)*(-0.5))*£(j); 
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end 
sum = sum + range(ring) *(y(9)*(-0.5))*£(9); 
temp = (step/3)*sum; 
end 
if z == 6, 
pr(ring,i) = temp * 2*(-3/2); 
elseif z < 13.7, 
pr(ring,i) = pr(ring,i)+ 2*temp*z* (-3/2) ; 


elseif z == 13.7, 
pr(ring,i) = pr(ring,i)+temp*z*(-3/2) ; 

end 

end 

pr(ring,i) = pr(ring,i)*.1/2; 

end 

% calculate p(3,i) for z = 13.7 to infinitive. 
if ring == 1, 


for i= 1:10, 
y = 1:.1:26.6; 
clear f; 
f = polyval(fit(i,:),y); 
sum = range(ring) * (y(1)*(-1/2)) * £(1); 
for j = 2:2:256, 
sum = sum + 4*range(ring) *(y(j)*(-0.5))*f(3); 
end 
for j = 3:2:255, 
sum = sum + 2*range(ring) *(y(j)*(-0.5))*£(4); 
end 
sum = sum + range(ring) * (y(257)*(-0.5)) *£(257) ; 
temp = (.1/3)*sum; 
temp = 2*(13.7%*(-1/2)) *temp; 
pr(ring,i) = pr(ring,i)+temp; 
end 
end 
% Throughput 
g = 0:0.05:4; 
for i = 1:81, 
s(ring,i) = 1/(2*range(ring) )*g(i)/((g(i)+1)%3); 
for n = 1:10, 
s(ring,i) = s(ring,i) + (n+1)*(n+2)/2*pr(ring,n) 
Z (SPEARS (SENG a ee CT) AB AOS Te) 
en 
end 
end 
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APPENDIX B. COMPUTATION OF CONDITIONAL CAPTURE PROBABILITY 


IN RING R, 


Let us consider a tagged packet in ring R; bounded by two 
concentric circle of radii r,; and r,,,, respectively. The 
conditional probability density function for packet generation 


in ring R, is given by 


£,(z) r,.<2r62 
fp (r|R,) = Fy (2341) =F, (r,) J itd (Bl) 
0 , otherwise 


The probability a packet is generated within a distance 


rsr,; is 


Zs 
2n G 
me a! n (R2 -R?) 


Solving (B2), we have 


F,(r;) (B3) 
2 ~Ry 
Similarly, 
2 
Tis 
Fe (2ieg) = =. (B4) 
2 ~ Ry 


Substituting (B3) and (B4) into (Bl), we obtain the pdf of 


packet generation in ring R, as 
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2zr 
fp (r|R;) a ria-ri 


0 , otherwise 


p AGS 2S Bey (BS) 


The corresponding conditional pdf for packet power X in 
ring R, can be obtained from (16) as 


3 
-(+) 
x 2 


2 (r?,,-277) 
0 , otherwise 


fy (x|R;) 0 fy, ST ST ia (B6) 


The pdf of the random power of n interfering packets 
anywhere in the network is obtained from (19). Therefore, the 


conditional pdf of signal to interference ratio in ring R, is 


£,(z|n,R,) = [yf,(yz|R,) £y(y|n) dy (B7) 
fe) 


Then the conditional capture probability in ring R, is 


Pr{capture|n, R;} = ff. (z|n, R,;) dz 
Yo 
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APPENDIX C. DERIVATION OF THROUGHPUT OF A RING IN A SYSTEM 


OF MULTIPLE RINGS. 


The derivation of the ring throughput for a multiple ring 
system is similar to the derivation in section II.B. 
The channel throughput S,, in packets per mean packet time 


becomes [5] 


+o; (x,u) a(x) a(u) dxdu (C1) 


s,-f 
ie) 


where a(x) and a(u) are exponential density functions with 


Om, 8 


mean x of the tagged packet length and the preceding packet 
length. The quantity ¢,(x,u)/2 is the conditional channel 


departure rate in ring R; [5]. 


SO, (x, u) = 59: (x+u) Pr(capture) (C2) 


where g, is a parameter in packets per second of a Poisson 
process of the channel rate in ring R;. 


Substituting (26) and (C2) into (Cl), we have 


= n! 
. Pr{capture|n} a(x) a(u) dxdu 


[$s (x+u) > e-gixeu) _[g(x+u)]” 2 
0 


Utilizing the binomial identity 


( yn = “jn k ,)n-k (C4) 
x+u Ye (al* u 
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and expanding (C3), we obtain 


25 gif ft xe"* e-% + ye" e- D> 2" pricapture|n) 

29 o 0 Bre (C5) 
z>. (i)*" u®-* a(x) a(u) dxdu 
k=0 


Further manipulating terms in the above equation, we get 


1° $9 ‘oe 2" pr {capture|n} 


1 + - n-ko- 
y(2 Ip te malnian, u®-Ke-%a(u) du (C6) 





k=0 
oe 


r u2-k+ - ko- 
ale xy? le-%a(u) du. ie. e-% a(x) dx 


We have assumed that the packet length is exponentially 


distributed and the average length is x. Hence 


a(x) = te-*/® (C7) 
x 


a(u) = +eV/F (C8) 
x 


Define the average traffic rate G of the system as 


G = gx (C9) 


We apply (C7), (C8), (C9) and the following equation [9]. 


|=? e- #* dx = nipo (ae) (C10) 
0 


to solve (C6) 


The first integral term can be evaluated as 
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7 7 -(g+%)x 
f x08) ema (x) dx =f x (k+l) 9 OI" Fax 
is) Qo 


= (k+1)! (gay 


(k+1)! ~(k+2) 
x7 (k+1) (G+1) 


Similarly, the second integral is found to be 


-(g+i)u 
fu e-ma(u) du 7 =f y (2-b) ee 
fo) 0 


= -(n-k+l1) 
= An-kK)! (9,2) - 
aoe x 
= (n- : -(n-k+1) 
hearer 


The third integral term is 


pac e-% a(u) du = (n-k+1) ! (G+1) > (2-k+2) 
; x (a k+d) 


The final integral term is reduced to 


f x"e"%a(x) dx = m8 (G+1) 74m) 
Q x” 





(C11) 


(C12) 


(C13) 


(C14) 


Combining all integral terms in (C11), (C12), (C13), and 


(C14) into (C6), we obtain 


n 


= n - (n+3) 
S; = uy £* Pri capture|n} (St ni(n+2) (C15) 


fer x (4) 


Applying (C9), the equality 


> n! (n+2) = n! (n+1) (n+2) 
k= 
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nd the fact that G, = g,;x, and further simplifying (C15), we 


obtain 


S; = > {a+i)(n*2) pr{capture|n} G,6"(G+1)- "> (C17) 


n=0 
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APPENDIX D. DERIVATION OF PROBABILITY OF CAPTURE. 


The probabilit: of capture in a ring Pr{capture}, per mean 
packet time can be found by 
Pr{captuwce} = 


ned 


a(x) a(u) dxdu 


1; Pr{capture|n} Pr{n} (D1) 
i?) 


oo, 8 


where a(x) and a(u) are exponential density functions with 
mean X of the tagged packet length and che preceding packet 
length. 


Substituting (26) into (D1), we have 


Pr{capture} = [ i> engixew Lox u) ]* (D2) 
0 ‘ono : 


. Pr{capture|n} a(x) a(u) dxdu 


Utilizing the binomial identity 


n 


(xeuyt = SP (a )xk ur (D3) 


and expanding (D2), we obtain 


oo Ld - 7 
Pr{capture} = f{x <~ pricapture|n} e7'(= 
0 o7=0 : (D4) 


3 (i)** u®-* a(x) a(u) dxdu 
k=0 


Further manipulating terms in the above equation, we get 
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= n 
Pr{capture} = > 2_ pricapture|n} 


n=Q 


(D5) 


ys ( 2) fxte-traton ax. f u®-ke-%a(u) du 


k=0 


We have assumed that the packet length is exponentially 


distributed and the average length is x. Hence 


a(x) = te-*/* (D6) 
x 


a(u) = tek (D7) 
x 


Define the average traffic rate G of the system as 


G = gx (D8) 
We apply (D6), (D7), (D8) and the following [9] 


ieey e- 8% dx = nipo (re) (D9) 
0 


to solve (D5) 


The first imtegral term can be evaluated as 


© st 
[ x*e- a(x) dx =f x*ke (9°) X ay 
0 


x 
1?) 


- (k+1) 
meget ip?) 
x Xx 


u 


k _! (G1) ~(**0) 
x* 


Similarly, the second integral is found to be 


/ ul?) e-%a(y) du = LAK)! (gy) - take (D11) 
3 a-k) 


ie) 











Combining all integral terms in (D10), (D11), into (D5), 
we get 


= n \ -(n+2) 
Pr{capture} = ye - pr{capture|n} $2*7) 
n=0 % 


-n 


Applying (D8) and the equality 


> nm! = n! (n+1) 
k=0 


and further simplifying (D12), we obtain 


Pricapture} = )* Pr{capture|n} (n+1)G7(G+1)~("*2)— (D14) 


n=0 
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